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Abstract 

Parton equilibration studies for a thermally equilibrated but chem- 
ically non-equilibrated quark-gluon plasma (QGP) , likely to be 
formed at the relativistic colliders at BNL and CERN is presented. 
Parton equilibration is studied enforcing baryon number conserva- 
tion. Process like quark - flavour interchanging is also taken into 
consideration. The degree of equilibration is studied comparatively 
for the various reactions / constraints that are being considered. 
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1 Introduction 



The partem plasma likely to be created in high energy heavy ion collisions is 
sure to be in a state of both thermal and chemical non-equilibria. Once thermal 
equilibrium is attained, the subsequent evolution takes place according to the 
laws of hydrodynamics. 

It is known that [ 1 ] the initial system produced at RHIC energies have 
a finite non- zero baryon density. The anti-particle to particle ratios at mid- 
rapidity for y / sjv : /v=130 GeV Au+Au collisions at the STAR collaboration, BNL 
report that there is a noticiable excess of baryons as compared to anti-baryons 
as reflected by the yields of ^ ( hovering between 0.6 and 0.8 ), A ( hovering 

between 0.7 and 0.8 ) [ 2 ] for varying transverse momentum, centrality and 
rapidity. In a subsequent reporting of the BRAHMS collaboration of BNL[ 3 ] 
it is seen that during Pb+Pb collisions at y/SNN= 200 GeV the observed 2 
ratio is a maximum of around 0.8 at zero rapidity and falls steadily for higher 
rapidity values. 

This excess of baryons ( and hence quarks ) over anti-baryons ( and hence 
anti-quarks ) clearly necessiates the inclusion of a chemical potential into the 
theoretical framework . Moreover, strong interaction conserves baryon number, 
which also must be incorporated in our scheme. The present piece of work is 
aimed at attaining such an objective. 



2 Brief Review of Previous Works 



Inclusion of chemical potential into the framework of Parton equilibration 
studies has been rather recent. Although in 1986 Matsui et.al reported [ 8 ] on 
strangeness production rates at nonzero chemical potential, it was not until the 
very end of the past century that a group from BARC,Mumbai, India included 
a chemical potential into chemical equilibration studies [ 4,5 ]. They considered 
a chemical potential (for massless quarks) that equalled the system temperature 
at all temperatures which was not realistic as it did not consider baryon number 
conservation. Furthermore, they did not distinguish between quarks and anti- 
quarks. The same criticism regarding chemical potential and baryon number 
conservation also applies to an attempt [ 6 ] (incomplete) by the present author 
during the preliminary stages of the current work. However, explicit masses 
and distinction between quarks and anti-quarks had been incorporated in that 
conference report. In 2004, He et.al. [ 7 ] undertook a much more complete 
study including baryon number conservation. They used an expansion of the 
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number densities over powers of chemical potential. We basically follow the 
same line with a few points in difference with the work by He et.al.The points 
of difference are: 

i) We explicitly introduce the strange quark mass 

ii) We treat quarks and antiquarks separately.lt can be noted here that al- 
though He et. al. [ 7 ] started off with quark-antiquark distinction, subsequently 
that distinction was put away. 

hi) We include results of full phase space calculations for both quark-antiquark 
pair production and quark flavour changing processes. 



3 The Chemical Potential 

At this point let us note the following points: 

1. The initial strangeness content of the QGP fireball is zero and moreover 
strangeness is very nearly conserved in strong interaction which clearly indicates 
that the strange quark chemical potential must be zero. 

2. If we consider a pair production process in QGP ( production of non- 
strange quark- antiquark pair) it can be argued that the net chemical potential 
on both sides of the reaction must be the same. Since gluons have no chemical 
potential clearly this would imply 

m = -\Xq 

...(1) 

3. We treat the up and down flavours in the same footing and treat them as 
massless. 

4. For a given initial baryon number density and initial values of the light 
quark & antiquark non-equilibrium fugacities ( which is defined as the ratio of 
the number density to the equilibrium number density) we can find the initial 
light quark chemical potential by method of iteration. 

5. The rate of change of chemical potential with time is dictated by the 
baryon number conservation equation. 
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4 Thermodynamic &: Statistical Mechanic Treatments 



4.1 Distribution functions 



For a thermally equilibrated but chemically non-equilibrated parton plasma, the 
distribution functions are given by the Juttner distribution functions [ 9 ] . In the 
present paper, initially we use the full distributions. In the latter stage we shall 
adopt certain approximations in order to achieve some mathematical simplicity 
as also to confront some mathematical unsolvability that we shall encounter. 
The full Juttner distribution functions are 



f — ^ 9 

Jt 

tri 

A 



e^-X g 



£ A l(l) /r,\ 

f =f-= 



4.2 Basic Quantities 



Using standard techniques, we can obtain the various thermodynamic quantities 
that are required to specify the state of the system. We obtain: 

Q = Q g +Zi= u ,dAQi + Qi) =Q g + 2*(Q q + Qq) + 2 Qs -( 3 ) 

Q being the number density, energy density or pressure in general and q denot- 
ing the massless quark flavours. 
This gives 

. . k/j.q , k/j,q 

16T 3 A™ 12T 3 . n X k „e~ + A£e ~ 

» = (— + — ^s-if- 1 ) <- — jH — )+ 

48T 4 A* 36T 4 , , A?Se"r £ + A^e r~ 
e = {-^5-=£i fcl + -y^(-l) k -H^ j?* }+ 

6m (L ^oo / n fc-i x fc/ g^(fcgg) . Kijkxs) 
— £ fc=1 (-l) A s { (fc ^ )2 + -^—^11 

and k 

16T 4 A^ 12T 4 t , A'fe^ 3 " + A^e t 2 " 

» = (— + — =£.(-!)'-'{-' ^ ) + 

...(4) 



4 



4.3 Baryon Number Conservation Equation 



The net baryon number of the system equals 

n B V = -{n q - n-q)V 
If we restrict our discussions to a unit volume then the baryon density equals 

— g iTi? - T&gJ — -^-2^=11.-1) { p } 

for given values of the parameters. Thus, if the initial baryon number density, 
temperature and non-equilibrium fugacities be known ( as from models like HI- 
JING or SSPC or likewise ) then we can iterate to get initial value of the quark 
chemical potential. 

The baryon number conservation equation gives 

(n B u") = fy* + ^ = ...(5) 
from which we can obtain expression for fi q , which gives the rate of change of 
the chemical potential with time. If we substitute for the number densities and 
solve for fi q we get 

M 9 =X q B 1 +X 1 B 2 +fB 3 + B 4 -..(6) 
where 

^i- x-^=i — p — \ e 



^ Z^fc=l P A 9 



e 



3 . , A„e~ + A£e r~ 

^3 = -^ 1 (-l) fc - 1 {^ ^ > T 

and 



rA 



b 4 =--i r sgL 1 (-i)*-H ^ J J } 

/, -I , V.' ' +->£e~ 



with 



a= ;r } 

We shall see the effect of having these nonzero values very shortly. 
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4.4 Quark &; Anti-Quark Number Density Evolution Equation 



Using standard techniques we can obtain the fermionic number density evolu- 
tion equations. The marked difference from the zero chemical potential case is 
that due to the nonzero values of the coefficients of the Baryon Number Con- 
servation equation, the massless quark and antiquark number density evolution 
equations both depend on either of the two growth rates. Let us put forward 
the two Number Density Evolution Equations. We include a new process called 
the quark flavour changing process (qfcp) into our framework. 



4.4.1 Massless Quark Number Density Evolution Equation 



If the Quark Number Density Evolution Equation be explicitly worked out it 
takes the form 

XjQi+Vfc + TQ 3 + Q 4 = -(7) 
where 

Qi = -^-e^C-i)*-^- 2 ^ ^} 

1W3 \k e ^B2 



, .AjT^-Tr 12T 3 , ,Aj;e^^ 

V°o / i \k — 1 9 \^oo / t\k— 1 1 T 

\ fc e T 19T 3 X k e^^ L 



Q 3 = 


36T 2 






IT 2 


and 




Qi = 


12T 3 






T7T 2 


where 





4.4.2 Massless Anti-Quark Number Density Evolution Equation 



For anti-quarks the equations are more or less identical as above with X q — > A, and 
fiq — > — /i g , except for the modifications in the contributions from Bi . If 

the Anti-Quark Number Density Evolution Equation be explicitly worked out 

it takes the form 
X q AQi+XqAQ 2 + TAQ 3 + AQ A = -(8) 



C 



where 



12T 3 , , Ke t - F j- 

AQi = ^^li-lf- 1 ^ 



k 2 



k[i q fe/i 



AQ 2 = ^S^t-^-H ^ - ^ } 
_ 36T 2 fc, ^*-^-^ 12T 3 Afej^f 

i kfj,q , kfj,q R 

12T 3 , . lL1 Ae-- 12T 3 _ , ., tJ Ae-T| 



^4 = ^E^-lj'-i-Lj ^S-Sg^-l)*- 1 f fc2 " ^ SaQ 

where 

S(lQ — i^Rgg — >g<j — ►fl'S/ (^ss — >gq — >ss) = 5Q 

4.4.3 Massive Strange Quark Number Density Evolution Equation 



For massive strange quark the evolution equation is given by 



X s S 1 +TS 2 + S 3 = ...(9) 
where 

1 — : — \ — 

s E^(-i)- 1 AS{^ + ^)> 

<->2 
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E *=l ( 1} A * (fen.) 



5.3 =- - 5"Qs ( He ) with 

T 

SQs = {(Rgg^sg — i? S s_+gg)/n s } + 2{(i?,jq_» S s ~~ Rs J — > <jg ) / n s } 



4.5 Energy-Momentum Conservation Equation 



The conservation of the energy-momentum tensor , which we term as Bjorken's 
equation for future references, gives 

^ + ^ = ...(10) 
Inserting the full expressions of energy density and pressure we obtain 
Tfe + Kjh + Va + X 7 / 9 + A./. + /10 = ...(H) 
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where 

192T 3 ~ X k g 36T 4 B 3 ^ A|e^-A|£^ 

fc=i 

144T 3 A fc e ^ + A | e -^ 72mg,^ 00 fc _ t fc if 2 (fcx,) 
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h-—^ k =i^) { p } + p 

_ 6rr4 fc _ x fc _ x 3X2^) ^(tes)-. 
Js — — 2~ n k=l\~ L ) K - A s \ — T, s 1 a \ / 

7T (KXs) (KX S ) 

64T 4 ~ x h n 48T 4 t , Ajfe^ + A£e~^ 
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fe=l 

2~^fe=ll--LJ * s \ Jl \ 1 7] >— > + 

T7T (KXs) (KXs) 

36T 3 . , \„e~ - A£e ~ 
^-Eg^-l *-!{-« rj 1 }B, 

5 Changing Distribution Functions 



In order to reduce computational complications to a low degree we make cer- 
tain approximations. We redefine the non-equilibrium distribution functions and 
shift from the Juttner distribution function to Modified Fermi-Dirac type for 
the quarks and antiquarks in confirmation with earlier works [4,5]. Needless to 
say, this would surely limit the applicability of the present study. However, it 
is believed that the basic nature of partonic fugacity variations would remain 
unaltered. Since we adopt certain approximations for reasons as mentioned ear- 
lier, we refrain from giving the explicit forms of the reaction rates presently and 
hold back the expressions thereof for now. 

As stated, we adopt the following approximations: 

We replace the original distribution functions by the following: 
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f - Ag 

■13 e e/T _ l 

_ \{q) e T ...(12) 

■/9W) - e e/T + l y ' 

f = f-= ^ 

J s 3 e e/T < i 

Of these, the last two are the Modified Fermi-Dirac ( MFD ) type Distribution 
functions as was used in [4,5] . 

As for the quarks, we take 

i) For s quark we take mass~ 150 MeV. 

ii) For lighter quarks since mass is much too less than the strange quark mass, 
we treat the mass of lighter flavours as zero. 

iii) Since the lighter quarks are treated to have the same mass, we neglect 
uu ^ dd type of reactions as they would be like elastic scattering processes 

leading to no net particle production. 

5.1 Basic Thermodynamic Quantities 



Using the techniques of statistical mechanics, we can obtain the number density, 
energy density and pressure as follows: 

7T Z7T 

~^ h i^r 

e = — X g + 4Q J {\e^' T + A ¥ e ^ T } + 

6"4 S ^ fc.j 3K 2 (kx s ) Ki(kx s ) 
tt 2 As fr{ [ ij { {kx s f + (kx s ) 1 

8tt 2 T\ 7ir 2 n f T 4 rx „ /T , _ u /r . 
P = ~^ A 9 + 120 {A ^ e + A « e * }+ 

6™; 



0s \ VV 



fe=l 

tuit/i x s = m s /T 

...(13) 

Using standard techniques we can obtain the quark and anti-quark number 
density evolution equations. As for the Juttner distribution case, the equations 
are coupled via the Baryon Number Conservation Equation coefficients Bi . 
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5.2 Baryon Number Conservation Equation 



Using the modified distribution functions we can obtain the coefficients . 
We obtain 



where 
Ai = Xie^ + A ? e~"T" 
A 2 = A^eT" - A-e^T" 

5.3 Quark and Anti-Quark Number Density Evolution Equations 



Just as in the case of the full Juttner distribution, here too the quark and anti- 
quark number density evolution equations are coupled via the coefficients Bi . 



5.3.1 Massless Quark Number Density Evolution Equation 



The equation is given by 



fi q = X q B 1 +X 1 jB 2 +TB 3 + B A 



...(14) 



where 




X g Q 1 +XgQ 2 + TQ 3 + Q i = 



...(15) 



where 
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SQ 



T 

{{Rgg- 



SQ with 

>qq — Rqq^ 



,)/n q }-{(R q 



ss Rss- 



l)/ n q} 



5.3.2 Massless Anti-Quark Number Density Evolution Equation 



The equation is given by 

XqQl+XqQ 2 + TQ 3 + Q 4 - ...(16) 

AQi = -f 

AQ A =- - SaQ with 

r T 

SqQ {(Rgg — >qq Rqq — >gg)l'^'~q\ {(Rqq — ts's R-s's — >q'q)/^''q\ 



5.3.3 Massive Strange Quark Number Density Evolution Equation 



For massive strange quark the evolution equation is given by 

\ S S 1 +TS 2 + S 3 = Q ...(17) 
where 

*-£ 

s j&^^M 

S 3 =- - SQs with 

T 

SQs = {(Rgg^ss — R S s~ tg g)/n s } + 2{(Rqq^> S s — R S s->qq)/n s } 
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5.4 Bjorken's Equation 



From Bjorken's equation we get on substituting for energy and momentum 
T/6 + A g / 7 + AJ 8 + A ¥ /9 + A s / s + / 10 = ...(18) 

where 



32tt 2 



h = -7^— A S T 3 - 



40 v 9 q T T 2 40 

72mo s , nk-if^j(^) , 5ifi(fca; 8 ) K a (kx s ) 



4 in/tj/ 



n 2 T " s ^ v L (fcr s ) 2 ' \2{kx s ) ' 12 

f7 ~l5 T 

6to^ y> fc _! 3Jf 2 (fc:r s ) ATi(fcx s ) 
/s tt 2 ^ l j 1 (fc^) 2 + (kx s ) ' 

/s = ^pi e T + (V T - V r )^7> 

7n 2 nf _M£ M« -H±^B2^ 

/o = ^p{e ^+(A 9 e^~A^e r ) — } 

_ 32tt 2 4 6mg s , fc _ x 4^(fcgj) ifiQx^ 

/l0 - 457 A * T + W AsSfc=l(_1) + Tfc^T 1 

+ ~120 ^(V t + A^e - ) + -^^(A^ - - A ? e - ) — 

Now, it is to be noted as the quark and anti-quark number density evolution 
equations and the baryon number conservation equation are mathematically 
very simply connected, so these two evolution equations linked by the four B 
coefficients along with the energy-momentum conservation equation with inputs 
from the gluon and strange quark number density evolution equations do not 
constitute a set of solvable equations as the system determinant vanishes. To 
make the system solvable, we expand the exponential term in the numerators 
of the first two B coefficients and thereby introduce some 'assymetry'. Strictly 
speaking, this however would mean not exact baryon conservation. 



5.5 The Partonic Reaction Rates 

5.5.1 The Gluon Multiplication Rate R 3 

The gluon multiplication rate has been calculated by Xiong et. al. [ 10 ]. The 
rate of the reaction gg — > (n — 2)g is given. By explicitly calculating the ma- 
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trix element [ 11] (summed over all the final states and averaged over all initial 
states) we can obtain the gluon multiplication rate. 

However, to avoid the huge calculations of evaluating 25 Feynman diagrams [ 11 
] involved, we fall back on previous results that earlier workers had used . We 
postulate that the gluon multiplication rate depends on the chemical potential 
via the Debye Screening mass. 

The Debye Screening mass suitable for a multicomponent chemically non-equilibrated 
parton plasma is given by [ 12 ] 



r> = \ I dkk[N c f g + J2i fi] -(19) 



m 

7T 



where the sum runs over all flavours i, while N c gives the number of colours. 
To accomodate for antiquarks and remembering that our number of flavours is 
3 and not 6, we propose the following modification: 

m » = %l dkk ^ + \ T,i=u, d Mi + ft)] -( 2 °) 

Using standard techniques , we can obtain the following results for the mean 
free path: 
\-i _ (a 2 dull _ QUg-KCil 

= %M ^%„4<1 + ^) ••■<"> 

which for zero chemical potential reduces to the well known result: 

x-i - 9 1 

X f -8 aiasT l + 8na s \ g /9 
Using methods of our previous chapter we get on integrating the modified dif- 
ferential cross section and recalling the definition of R 3 = -a^Ug the required 
rate as 

f = 2 -^X}ngI(Xg) -.(22) 
where 

s 

J(A-) = / dx 4m c dz- -[ V + 

9 Ji Jo (1 + z) 2 l Xy /\x + (1 + z)x D } 2 - Axzxd 

1 coshr X \fx 



sX j + x(l + z)y D ] 2 - 4xzy D 

with 

xd = m 2 D Xf 

m D 
VD = 

s 

5.5.2 The Quark Anti-Quark pair production rate 

We have parton production rates in the RHS of the number density evolution 
equations. Let us evaluate the quark-antiquark pair production reaction rate 
R 2g . We have [ 8 ] 

R2 9 = R% n Rf£ s -..(23) 
where 
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9 J (2tt) 3 2E 1 J {2tt) 6 2E 2 J (2ir) 6 2E 3 J (27r) 6 2E 4 

6 4 (Pi+P2-P3- Pi mM 99 ^l\ 2 f 9 (Pi)fM)(l ~ /,(Ps))(l ~ f'M) -( 23a ) 
and 

d 3 Pl f d 3 p 2 f d 3 p 3 f d 3 pj 



RL = —4^- —4^- 



l ° SS ' (2tt) 3 2E 1 J (2tt) 3 2E 2 J (2ir) 3 2E 3 J (2tt) 3 2E 4 

s 4 (p 1+ p 2 -p 3 ~ Pi mM gg ^ a \\i + / s (pO)(i + f g ( P 2))f q (p 3 )fM -( 23b ) 

Following [ 8 ], we can say that there are three topologically distinct Feynmann 

diagrams that contribute towards the quark-antiquark pair production process. 

Evaluating them performing traces and finally adding them up we can find the 

net squared matrix element. We basically follow the lines of [ 8 ]. Transforming 

variables as 
q = pi + p 2 

P = jj(Pi -Pz) 

q =P3+ P4 

P = -Pi) 

with restrictions 

qo > 2m,i 

« = «g-|qf >4mf 

and transforming the three dimensional integrals to four dimensional integrals 
using 



d 3 Pi 



2\ 



2E- I d4piS (P 2 ~ m *> 
with the new set of variables 

qo = —Tlnv + 2rrii 

= {ql-im 2 i)iu 

q , 4m? . i 
Po = |(l % -)*x 

we arrive at the rate 

^ /"I /"I /"I „2 Arm 2 



or f f f f u l Am 2 i 3 

R2g = 7T^ T du dv dx dy—(l -) * (gg - 4m- ) 2 f Quarks fphasel -(24) 

Jo JO Jo JO u s 

where 

/Quarfcs = /<,("£ + Po)fg(-f ~ Po)(l ~ fqif + Po))( 1 ~ ~ Po))' 



f +P0 ))(1 + / 9 (f -po))/,(f +Po)/?(f 



and ^ ^ 

/ pftasel = A + B[— + —}+ C[j^ + 

with 
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. oM M 4m? .(1 -x 2 )(l -y 2 ) 2 2ll 34 n m 2 

A = 3 1-1 — + x 2 y 2 }} 24— 

s 2 3 s 

3 s s z 

128 mf 

tf± = [1-[1 - ^][(1 - ^ 2 - 2/ 2 ) ± 2[1 - ^Msy]* 



5.5.3 The Quark flavour changing rate 

For the quark flavour changing process we have [ 8 ] 

r w = r^-rZ s -( 25 ) 

where 



' J (2it) 3 2E 1 J (2ir) 3 2E 2 J (2tt) 3 2E 3 J (2ir) 3 2E 4 

2 



s(pi + P 2 -pa -p 4 )s|m s3 ^| f q ( Pl )M P2 )(i - f s ( P3 ))(i - fM) -( 25a ) 

and 

J (2it) 3 2E 1 J (2tt) 6 2E 2 J (2ir) 6 2E 3 J (2n) 6 2E 4 

SHP1+P2-P3 -Pi)nM SJ ^ q /(l - /,(pi))(l - U(p2))fs(p S )fM -( 25b ) 
Following [ 8 ], we can say that there is only one type of topologically distinct 
Feynmann diagram that contributes towards the quark flavour changing process 
or the strange quark pair production process. Evaluating it, performing trace 
calculations we can find the squared matrix element. We basically follow the 
lines of earlier works. Transforming variables as in the case before and perform- 
ing identical operations we can get the rate as 



^ Jo Jo Jo Jo v 



) 2 (<Zo — 4m?) 2 f Strange fphase2 



with 

fQuarks = /q(y + Po)fq (y ~ Po) i 1 ~ /«(y +P ))( 1 ~ /s(y-Po))- 

(1 + /.(f +Po))(l + /.(f -Po))/,(f +Po)/«(y -Po) 

4m- lr (l-i 2 )(l-w 2 ) , 4m 2 
= [1+[1 L ]P 77 — + aV]] + — 
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6 Results 



The initial Conditions from the Self-Screened-Parton-Cascade (SSPC) Model 
are as follows: 





Ti(fm) 


T(GeV) 




Xq 




A s 


RHIC 


0.25 


0.668 


0.34 


0.064 




0.032 


LHC 


0.25 


1.02 


0.43 


0.086 




0.043 



In absence of known value of initial baryon number density at the point when 
thermal equilibrium is attained (after which point the system evolves according 
to the laws of hydrodynamics) we perform a study with varying initial condi- 
tions. We vary the initial baryon number density between 0.11 and 0.21/fm"3 
while we change the initial quark to antiquark fugacity ratio between 1.1 and 
1.9. Needless to say, most set of values would be of academic interest only and 
not physically realisable. But, this could lead to a better understanding, it is 
believed. 

We study comparatively the outputs with given initial baryon number den- 
sity and initial light quark to anti-quark fugacity ratio. We arrive at the follow- 
ing conclusions : 

1. The nature of variations of the physical quantities are more or less in 
the line of earlier works. The variations of Temperature, Chemical Potential 
and Non equilibrium Fugacities for RHIC and LHC initial conditions are shown 
in figures 1 and 2 respectively. The decaying curves give temperature varia- 
tions while the positive rising curves are sequentially (from top) for gluon, light 
quark, light anti-quark and strange quark fugacity variations. Plots for both 
quark flavour changing process included and excluded cases are shown. For 
identification please see point 4 below. As representative values we take baryon 
number density 0.15 and initial fugacity ratio 1.5. 

We observe that 

i) As expected, the temperature falls with time while the non-equilibrium fu- 
gacities increase. 

ii) Contrary to the Juttner case, the chemical potential remains negative all 
along and as expected approaches zero as the system equilibrates. 

iii) The QGP,as expected remains to be gluon dominated. 

iv) For a system at higher chemical potential, the temperature falls at a slower 
rate signifying lesser amount of energy expenditure to create partons which 
shows up in the slower rise of all partons except for the light anti-quark, which 
shows a larger growth rate due to the presence of the exponentiated chemical 
potential. 

2. For a given initial ratio, except for the light quarks and anti-quarks, the 
output does not depend much on the initial baryon number density. For the 
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light quarks and antiquarks this variation is due to the presence of the expo- 
nentiated non-zero chemical potentials. 

3. For a given initial baryon number density we can recast the equation for 
baryon number density in the form 

.-(27) 

where C is a constant for a given temperature, light quark fugacity and given 
baryon number density. Here D is the light quark to antiquark initial fugacity 
ratio. For a fixed light quark fugacity and temperature, as D falls clearly the 
RHS of the above equation increases which indicates a rise in the chemical po- 
tential. Again for a system of higher chemical potential, the temperature has 
to drop at a slower rate due to the constraint imposed. Hence we observe that 

i) The temperature falls at a slower rate for a smaller value of the light quark 
to antiquark initial fugacity ratio. 

ii) As the temperature falls at a slower rate, it would imply a lesser expenditure 
of energy to produce particles in general . This would show up in the slower 
rise of all fugacity values except the light antiquark, which would show a higher 
growth rate. This is due to the exponentiated chemical potential part. 

4. For inclusion of the quark-flavour-changing process we observe the fol- 
lowing: 

Due to the additional production of s-quarks from massless quarks via qfcp, 
we see an additional increase in fugacity of the strange quark while, the rate of 
equilibration falls for non-strange fermions . 

5. For some initial baryon number density and known values of non-equilibrium 
fugacities and temperature switching distribution functions would mean a shift 
in the value of the light quark chemical potential, which is found by iteration. 
Since everything else remains the same, to fit the baryon number density, clearly 
the chemical potential is expected to 'pay the price'. This is exactly what hap- 
pens. For the MFD case, the chemical potential assumes a value much smaller 
(-ve and also magnitudinally an order smaller) than that for the Juttner case 
(in which case it is +ve) . In the result for the number density using the Juttner 

\ p — X* — I \ p — T — 

distributions we had an expression like ^oo ( l)*~ 1 f g " r/ ? i . 

If we expand the exponential with an objective of suitable truncation, keep- 
ing a 'k' factor in the numerator of the exponent with a much larger value of 
the chemical potential, it would make the mathematics still more complex as 
it would be required to consider many more number of terms for the same ac- 
curacy. Therefore, we had to adopt the MFD type of distribution functions ( 
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with k=l in the exponent ), where, as required, we can easily use a truncated 
expansion of the exponential term. 
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Figure!: RHIC Case 
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